ABSTRACT
Objective To screen the potential key genes of osteosarcoma by bioinformatics methods and analyze their immune infiltration patterns. Methods The gene expression profiles GSE16088 and GSE12865 associated with osteosarcoma were obtained from the Gene Expression Omnibus(GEO),and the differentially expressed genes(DEGs)related to osteosarcoma were screened by bioinformatics tools.Gene Ontology(GO)annotation,Kyoto Encyclopedia of Genes and Genomes(KEGG)pathway enrichment,and analysis of immune cell infiltration were then carried out for the DEGs.The potential Hub genes of osteosarcoma were identified by protein-protein interaction network,and the expression of Hub genes in osteosarcoma and normal tissue samples was verified via the Cancer Genome Atlas(TCGA). Results A total of 108 DEGs were screened out.GO annotation and KEGG pathway enrichment revealed that the DEGs were mainly involved in integrin binding,extracellular matrix (ECM) structural components,ECM receptor interactions,and phosphatidylinositol 3-kinase/protein kinase B(PI3K/Akt)signaling pathway.Macrophages were the predominant infiltrating immune cells in osteosarcoma.Secreted phosphoprotein 1(SPP1),matrix metallopeptidase 2(MMP2),lysyl oxidase(LOX),collagen type V alpha(II)chain(COL5A2),and melanoma cell adhesion molecule(MCAM)presented differential expression between osteosarcoma and normal tissue samples(all P<0.05). Conclusions SPP1,MMP2,LOX,COL5A2,and MCAM are all up-regulated in osteosarcoma,which may serve as potential biomarkers of osteosarcoma.Macrophages are the key infiltrating immune cells in osteosarcoma,which may provide new perspectives for the treatment of osteosarcoma.